Six-loop Konishi anomalous dimension from the Y-system 
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We compute the Konishi anomalous dimension perturbatively up to six loop using the finite set 
of functional equations derived recently in [1] . The recursive procedure can be in principle extended 
to higher loops, the only obstacle being the complexity of the computation. 



I. INTRODUCTION 

Using integrability in conjunction with the AdS/CFT 
correspondence led to some of the main achievements in 
theoretical high energy physics of the last decade. One 
of them is determining the spectrum of the anomalous 
dimensions of the Af = 4 SUSY gauge theory, or equiva- 
lently finding the spectrum of string theory in AdS§ x S 
background. For a recent review on the subject see [2]. 
For operators with large charges, the anomalous dimen- 
sion is determined by a set of Bethe Ansatz equations, [2], 
while for operators with small charges, the information 
about the spectrum is encapsulated into an infinite set of 
functional equations derivable from the thermodynamic 
Bethe Ansatz (TBA) equations [3] and known under the 
name of Y-system. For the AdS / CFT integrable system, 
the Y-system was conjectured in [4] and derived from the 
the TBA equations in [5] . The anomalous dimension of 
the Konishi operator, the shortest operator not protected 
by supersymmetry, is used as a testing ground for these 
equations, both for analytical and numerical computa- 
tions. In particular, in perturbation four [4] and five 
loop [6, 7] corrections to the Bethe Ansatz were com- 
puted using the Y-system, and they were found to co- 
incide both with the corresponding Liischer corrections 
[8, 9] and with the four-[10] and five-loop [11] pcrturba- 
tive gauge theory computations. At strong coupling, the 
result of the extrapolation to short, Konishi-like opera- 
tors [12] agrees both with the numerical results [13, 14] 
and with the string predictions [15]. 

Recently [1, 16], the Y-system was reformulated in 
terms of a finite closed set of functional equations. In the 
present work, we set up a recursive procedure to solve in 
perturbation the equations of [1] and we perform the ex- 
plicit computation up to six loops. The procedure can 
be in principle continued to higher loops, in particular to 
the double- wrapping order at eight loops. Double wrap- 
ping was already attained for the ground state energy in 
the twisted AdS/CFT [17], however our case is signifi- 
cantly more complicated because of necessity to account 
the displacement of the Bethe roots. 



II. THE SET OF FUNCTIONAL EQUATIONS 

The wrapping corrections are encoded in a finite set 
of functional equations [1] which were derived from the 



AdS/CFT Y-system by solving the Hirota equation in 
the infinite wings of the T-hook defined in [4]. The dif- 
ferent solutions are then glued together using analyticity 
constraints which reflect the physical properties the Y- 
systcm has to satisfy. We give below a brief summary of 
this set of functional equations. 

a. Input data. The functional equations depend on 
a specific operator through the position of the Bethe roots 
uj, via the objects B {±) = ]jf =1 (± - xJ),R (±) = 

Iljli \[~^ ( x ~ an d the Baxter polynomial Q ± — 

V X 3 

( — l) M f?( ± )i?( ± ). Here we use the conventions /I a l = 
f(u+ia/2), /± = /I* 1 !, andi(fi) = u/2g+i v / l - u 2 /Ag 2 
is the Jukovsky variable in the so-called mirror regime, 
with a cut on Z = (— oo, — 2g]U[2g,oo). For the Jukovsky 
variable in the physical regime, with a branch cut on 
the interval Z = \—2g,2g), we use the notation x(u) — 
u/2g + yJu 2 /Ag 2 — 1. By convention, we denote with 
a hat the quantities which depend on x(u), if we want 
to emphasize the position of the branch cut. For the 
Konishi operator there are two magnons, M = 2, with 
Mi = -t*2 = l/\/l2 + 0(g 2 ) and Q(u) = u 2 - u\. 

b. Parameterization of the T- and q-functions. 
When considering a state from the sl(2) sector, as it is 
the case for Konishi, the Y-system can be determined 
from only three functions: two real-valued densities with 
a finite support on Z, p and p 2 , and a complex- valued 
function U(u) analytic in the upper half-plane. A rela- 
tively simple formulation of the set of functional equa- 
tions can be obtained using T-functions in two different 
gauges denoted T a , s and £^ a ,s- The first gauge gives a 
simple solution of the Hirota equation in the right band 
s > a 

To, s = 1 , T hs = s + K s * P , T 2 , s = 7g' ] 71,7 s1 , (1) 
with (JC * f)(u) = ^ J^dv f(v)/(v -u) and JC S = 
£M _ &~ s l The second gauge gives a solution of the 
Hirota equation in the upper band of the T-hook, a > \s\, 

X,2 - q l ^ a] q l i a] , = (U^m- a ^) 2 X,i , (2) 

= q^ a] q [ 2 - a] + & a] q [ r ] + qt ] fc a] + Q^t^ > 

• J afi — <7i2 <Zl2 ?34 *?34 

_ [+a] -[-a] _ [+a] -[-a] _ [+a] -[-a] [+a] -[-a] 
9l4 9l4 ?23 ?23 9l3 ^24 ?24 ?13 ' 

We will also use the combinations s defined by (2) for 
a < \s\. The q- functions related among themselves by 
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9i = 1 , q 2 = -iu + JC * p 2 - JC*W pv 
q 12 = (u - ui - a)(u + u x + a) = Q + Sq 



{913,923 = 9l4,924> 



934912 
939l2 = 
W a 



qu 

= 913924 

92913 

[+"] -[- 
% 94 



E 



fc=0 
914923 , 



t/ 2 {l,92,9 2 2 } 
9j t 29i" 2 



[2fe+l] 



the Pliicker relations [1] are determined by p 2 and U as 
follows 

(3a) 
(3b) 

(3c) 

(3d) 
(3e) 

W = W . (3f) 



9l4: 



94912 



9292^3 



924 : 



[+»] -[" 

94 9 3 



Above, IU pv = i(T^[+°l + The definition of 

/£>2(w) is slightly different from the one in [1], such that 
here P2(u) is of the square root type, in the sense that 
P2 (it)/ \/4<7 2 — "" 2 is analytic in the vicinity of the real 
axis. 

c. Auxiliary integral equations. In the intermediate 
steps of the computations one needs to compute 3 quan- 
tities, Y\ t \, Y 2 .2, h. Y1.1, F 2 ,2 are determined from the fol- 
lowing relations, considered at Im(u) > 0: 



In 



Y i, 1-^0,0 
1-2,2^1,0 



1 #1,0 



x — a; 



^1,1^2,2 ■% 



UK * Im In 



= 2/C*Re 



0,0 



^2 
^1,1 

^1,0 



1 [1-0] 



a; — x 



In 



^0,0 



, (4a) 

[-0] 
.(4b) 



The function h is found from equations (5.38) and (6.14) 
in [1]. Only its large- volume asymptotic solution (11) is 
needed for the six-loop computation. 

d. Equations for p,p2,U. After Yi,i,l2,2j^ arc 
found, the set of functional equations can be closed by 
finding p, p 2 from 

1 + ^2,2 _ ^1,2^1,0 1 + Y 2t2 



-f[l+0]^[-l-0] 
'1,1 '1.1 



1 + Y,-! ^ M ' 1 + Yi.i T+.T^ 

equations valid for ueZ, and determining U from 

2 r T [21 

yi,i^,2^o:o n 



(5) 



hum 



Yi,i%: 

Y 2 ,2^i,0 



■^,i7Ti 



(6) 



.7"i,2^7i 

for lm(u) > 0. The solution of (5) has a one-parameter 
ambiguity fixed as explained below. 

e. Supplementary constraints. Let us now empha- 
size the role of the Bethe roots. Bcthc roots appear as 
the zeroes of the following functions: 

^10 ( u j) — , fixes the value of a , (7a) 

3\ t \{uj) — , fixes the ambiguity in (5). (7b) 

Using them and the Hirota equation ^i + ^i~ = 

^0,0^2,0 + {U + U~) 2 ^i, we conclude that ^,0=^2,0 
should have a double zero at Bethe root which appears 
to be a double zero of $o,o- The overall normalization of 
U is not fixed by (6) and should be defined from 

^1,1^2,2 



UU 



ore 
^0,1 



U € Z. 



(8) 



/. Exact Bethe equations. The set of equations 
above has a solution for a range of values of u\. To get 
the correct answer for the energy, one has to insert the 
value of ui which is fixed by the exact Bethe equation, 
which can be for example written [1] in the form 



hr 

h+ 



1 2,2 '1, 



r i-2] 
^ = 1 



y- n— -7-1+2] 

J 2,2 '1,2 lii 



at u = v,i 



(9) 



Asymptotic solution. The equations listed above de- 
pend on the parameter L which sets the large u behavior 
of h. In the large volume limit L — > 00 the Y-system, 
hence the functional equations above, can be solved ex- 
plicitly [4]. All the functions 9^, with the exception of 
912, and 93 and 94 are suppressed at least by a factor 
x~ L with respect to 91, q 2 and q 12 , so they are zero in 
the asymptotic expressions of 5F a ,s- 92 and q i2 are in this 
limit given by (92) as = -iu + JC * p 2 , (9i2)as = Q ■ The 
asymptotic values of p and p 2 are 



VW 2 



E as + 2 



(p 2 )as= -4 



VW 2 



E,, s — 2 



(10) 



where E as is defined in (21). One can check that both 
the equations (4) and (5) are satisfied by (Yi iY 2 2) = 



asymptotic values of h and U are given by 



X 2 



#+1 



Ah 




2 Jx(x,x j ) 



n 



=i e 



j=i 



where for x(x, y) one can use the BES perturbative 
expansion [18]. The shift operator D is defined such that 
D/ = f D = f + - The overall normalizatio n A^ is irrel- 
evant, cf. (9) and (6), whereas Ajj = \yE as {E as — 2) 

x exp(6 7ff 2 (l - 4 3 2 + 28 5 4 ) -I8C3/ +24. 9 6 (3C 3 + 5Cs)) 
is fixed by (8) with 7 the Euler constant. The term 
containing 7 depends on the regularization scheme for 
diverging sum. We we use the regularization prescription 

T^i = -i^{-iu). 

Although for Konishi-likc operators L is not large, the 
asymptotic solution is still valid up to at least L loops 
(four loops for Konishi), so the weak coupling is effec- 
tively a large volume limit. We use as a constraint that 
the full solution should reduce to the asymptotic one at 
weak coupling. In the following, we will continue to call 
'asymptotic' the above quantities evaluated at the exact 
position of the Bethe roots Uj, determined by equations 
(9). These quantities will therefore incorporate part of 
the wrapping corrections via the corrections to the Bethe 
roots. 



III. WEAK COUPLING EXPANSION 

The strategy of solving perturbatively the functional 
equations is to subtract from the exact equations the 



3 



asymptotic ones and use that the deviations of the T 
functions from the asymptotic values are small. The re- 
sulting equations depend on the ratio of the exact and 
the asymptotic values, so we denote (T) r = T j (T) as . The 
following quantities enter the functional equations 



H = ln\ 




r = ln 



Ti- 



lii I 



.(12) 



All these quantities are small in perturbation, and this 
will allow performing the expansion of the functional 
equations. For example, up to seven loops, one has 



H 



Sqi2 Sqi2 (?24 + 2< ?14<?2 + qi^2)U 



Or 



(Q + ) 2 



(13) 



The quantities U and q 2 arc determined asymptotically 
from (3a) and (11), and at the leading order they are 
equal to 



TP 



2g 4 



i 1 

3 u 



(14) 



As for the q- functions 913,^23 = <Zi4,<724, they are given 
by the sums in (3c) and in perturbation they have an 
array of equidistant poles which will be responsible for 
the appearance of the zeta functions in the final answer. 
For example, at the leading order, we have 

q 13 =18 g 4 (-iu + Q V> (1) H" + l/2)) , (15) 

where ip' n ) denotes the n-th derivative of the digamma 
function. The functions q 3 and 54 have a similar struc- 
ture. (734 is given by a double sum and it contributes 
only at eight loops and higher. A priori, the infinite sum 
in (3c) creates also infinitely many poles at positions of 
the shifted Bethe roots. However, because of the equal- 
ity {U+/U-) 2 = -Q[+ 2 1/Q h21 which holds at the Bethe 
root at first four nontrivial orders and which is just the 
asymptotic Bethe equation, the poles in 913 cancel out 
pairwise (except for the first one, which is cancelled by 
an overall factor Q). The cancellation mechanism still 
holds for qi4 and 924 because (q^ + q^ jas = (=^i,i)as and 
&i,i(ui) = 0. At least up to seven loops the g-functions 
are given by linear combination of the multiple Hurwitz 
zeta functions which we define by 



Vai, 



00 / 1 \ ai 

»«(") =J2 ( ^fJk ) W- ■«»(« + *( fc + !))> 



fc=0 



Va(u) 



(o-l)! 



^-^(-iu), a>l, (16) 



with coefficients which are rational functions of u . 

The algorithm for computing the perturbative compu- 
tation is summarized in figure 1. The interior loop deter- 
mines the densities p and P2, or rather their variations 
with respect to the asymptotic values, Sp and 5p2, from 
equations (4), analytically continued to u € Z, and from 
(5). In order to preserve the square root structure of the 
two densities, we need to make a rescaling z = u/2g. 



Run once for 
four orders in g 2 




(yi,i"y 3 ,j) 



Run for each order in g 2 



FIG. 1. Structure of the perturbative computation. 



When performing various integrations, we encounter two 
different situations. In the first, the integration is (par- 
allel to) the real axis and we deform the contour to pass 
just below it and avoid the possible singularities on the 
real axis. This allows us to perform uniform expansion 
in g. In the second situation the Cauchy kernel K, acts 
on a function with a finite support, like bp and 5p2, then 
it can be expanded in terms of moments: 



n>0 V J 



except for s = 0, where /C[ ±0 1 * p — ±p/2 + |C * p , the 
slash meaning principal value. The equations (4) and (5) 
finally reduce to the linear system 



2 12 
1 3 



g 2 v 2 



9 1 Co 



(18) 



where C\ and C2 are Taylor series in g 2 and z is known 
from the previous orders of perturbative expansion (at 
the leading order in g they are constants) and 

Sp f(l-2z 2 )Sp 2 



^/Y^z- 



-r, V2 



Aiz fc*5p 2 .(19) 



By definition of p (and P2), we know that p(2gz) / \/l — z 2 
is analytic in the vicinity of the real axis and hence can be 
Taylor-expanded. At a given order in <?, this means that 
Sp and Sp2 are polynomial in z times yl — z 2 . Using this 
information, the equation (18) gives, at leading order, 

Sp 2 =g 7 {A + Bz 2 )^\-z 2 , 5p = g 9 C\f\ - z 2 , (20) 

with &A = Ci - 2C7 2 , 2C = -Ci + 4C 2 . The value of B 
is unconstrained by (18), and it is fixed by (7b). 

The result for the energy. The energy can be 
computed [1] from the behavior at large u of \nYi xY 2 _2, 
lnYi^Y^ — iE/u. Isolating the asymptotic and the 
wrapping part in the above expression, E = E as + £ , W rap, 
one obtains 



E as = 2-8<7lm 



1 



EL 



-H(u)du 



1 



hit 



(21) 



The Asymptotic Bethe Ansatz [18] predicts up to 6 loops: 
£bae = E &s {ui -> ui,bab) = 2 + 12 5 2 - 48. 9 4 + 336 5 6 
-(2820 + 288Cs)5 8 + (26508 + 4320 C 3 + 2880 Cs) 9 W 
-(269148 + 55296 ( 3 + 44064 C 5 + 30240 (7) g 12 ■ (22) 
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At finite volume (L = 2 for the Konishi operator) , the en- 
ergy receives corrections both from E as through the cor- 
rection of the position of Bethe roots U\ , and from E wrap . 
The corrections to the Bethe equations leading to dis- 
placement of the Bethe roots start at 5 loops, where they 
are due to correction of only. At 6 loops one should 
take into account the first corrections to p 2 , whereas ef- 
fects from correction to p and h are delayed at least up 
to 7 loops. -Ewrap is non-zero starting from 4 loops. The 
single- wrapping corrections, up to seven loop, follow from 
computations within interior loop in figure 1, whereas ex- 
terior loop has to be run only once, to find the explicit 
analytic expression (11). For H one can use the approx- 
imation (13). The double wrapping effects, in particular 
correction to (?7) as , are important starting from 8 loops. 

We have performed the explicit perturbative expan- 
sions discussed in previous sections in order to compute 
5E up to six loops, i.e. up to g 12 term. Intermediate ex- 
pressions are too bulky to be presented here. We summa- 
rized them in the Mathematica notebook file [19]. These 
expressions contain ^-functions (16) and their residues at 
the Bethe root. However, the final expression is signifi- 
cantly simpler and is given in terms of zeta-functions: 

SE 4Sz5 loop = (324 + 864 Cs - 1440 C 5 ).g 8 + (23) 
(-11340 + 2592 C 3 - 11520 Cs - 5184 Cf + 30240 Ct)5 10 , 

$E 6 i 00 p = (261468 - 207360C 3 - 20736C| + 156384C 5 
+155520CsC5 + 105840Ct - 489888Cg)g 12 • (24) 

At four and five loops we reproduced the already known 
answers [8, 9]. We were informed that Z. Bajnok and 



R. Janik have obtained the same result as us at six loops 
[20] using Liischer corrections. Also, fitting the known 
numerical results [13, 14, 21] with a diagonal Pade ap- 
proximant, we were able to fix the 6- loop energy with 5% 
confidence. Our analytic result is compatible with this 
numerical estimation. 



IV. CONCLUSION 

We have computed the wrapping corrections of the 
Konishi operator up to 6-loop order using the functional 
equations proposed in [1]. We have adjusted the struc- 
ture of the functional equations for a systematic pertur- 
bative expansion and we hope to be able to apply our 
methods to reach double-wrapping orders for the Kon- 
ishi states. An interesting question to explore is whether 
the cancellation of poles at Bethe roots observed when 
computing q t j is present at any order and if it can be 
related to the exact Bethe equation. Another question 
is what type of functions can appear in the final answer. 
So far the expression for the energy is reducing to Eulcr- 
Zagier sums and we believe that it will be always so. 
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